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Abstract. A procedure is described of construction of velocity fields and velocity dispersion of 
the stellar component in the central regions of galaxies on the basis of data obtained with the 
integral field spectrograph MPFS. Cross-correlation and Fourier quotients techniques adapted for 
2D spectra reduction are used. A detailed discussion is given of the problems of providing for 
the instrumental profile in measuring the velocity dispersion. Using the existing version of the 
spectrograph, one can measure radial velocities and the dispersion of velocities to an accuracy of 
5-10 km/s and 10-15 km/s, respectively. As an example the "double-barred" galaxy NGC2950 is 
considered. 
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1. Introduction 

The present paper examines the problem of deter- 
mination of two basic observational parameters that 
characterize the collective motion of stars in galax- 
ies — radial velocity v and velocity dispersion a. 
The task has a half-a-century history, beginning with 
the pioneering observations by Minkowski (1954) to 
the latest papers concerning the relation between the 
mass of a central black hole and the velocity disper- 
sion (Gebhardt et al., 2000). The procedure of mea- 
suring V and a with the application of the technique 
of Fourier transform goes back to the classical works 
of Tonry and Davis (1979, hereafter TD) and Sargent 
et al. (1977, hereafter SSBS). 

The methods of two-dimensional (2D) spec- 
troscopy that have been advanced intensively over 
the last decade have made it possible to obtain two- 
dimensional mapping of velocities and velocity dis- 
persions. In the case of disk galaxies this enables one 
to study non-circular motions in the region of the tri- 
axial potential (bar) or near the active nucleus (see 
the survey of Arribas and Mediavilla, 2000). The mul- 
tipupil spectrograph (MPFS) has been in operation 
at SAO RAS since 1990. It has been used to con- 
duct similar observations (see, for instance, Afanasiev 
et al., 1996; Sil'chenko, 1998). In these papers, how- 
ever, field velocities of stars and gas alone have been 
discussed. At the same time, the distribution of ve- 
locity dispersion gives information about the non- 



axisymmetrical shape of the gravitational potential. 
This is one of the key parameters in producing numer- 
ical dynamical models of galaxies (Khoperskov and 
Chulanova, 2001). The main methodological problem 
in measuring velocity dispersions is taking correct ac- 
count of the variations of the instrumental profile of 
the spectrograph. A new version of the MPFS (in op- 
eration at BTA since 1998) constructs reliably the 
distributions of velocity dispersions in circumnuclear 
regions of galaxies, which is shown in the present pa- 
per. 

The paper is organized as follows. The theory of 
methods of estimating v and a is considered in Sec- 
tion 2. The problems of provision for the instrumen- 
tal profile are discussed in Section 3, Section 4 gives 
a brief description of algorithms of analysis of MPFS 
spectra and errors that arise, the results of the galaxy 
NGC 2950 are presented in Section 5, conclusive re- 
marks are collected in Section 6. 

2. Methods of velocity and velocity dis- 
persion measurement 

2.1. General remarks 

After a preliminary reduction, the MPFS observa- 
tional data can be presented as two-dimensional spec- 
tra. The spectra from different spatial elements are 
brought to a common wavelength scale and arranged 
in the frame sequentially one above another. An ex- 
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Figure 1: Example of MPFS spectra: the 2D spectrum 
of the galaxy NGC2950 (bottom); the spectrum cor- 
responding to the galaxy centre (top). 

ample of such spectra is displayed in Fig.^. The spa- 
tial elements are constructed as small lenses and 
grouped into an array of 16 x 15 pixels (micropupils). 
This corresponds to a field of view of 16" x 15" on 
the sky. 

Let us illustrate the idea of the methods of mea- 
suring V and <j by an individual spectrum of one 
spatial element. Let G{\) and T(A) are the ob- 
served spectra of the galaxy and a template star 
with the low-frequency component (continuum) sub- 
tracted. Typical examples of such spectra are shown 
in Fig.||. The spectra are discretely sampled in wave- 
length A. 

The observed spectra need to be converted to a 
logarithmic scale in which the number of the spectral 
bin n is related to the wavelength A by the expression 



n ~ alog(A) + h, 

where 

a = log"^(l + Aw/c), 
b = -a log(Ao), 



(1) 



(2) 



here Aw is the step in velocity expressed in km/s, c is 
the speed of light. The Doppler shift leads to uniform 
shift of the spectrum along the n axis, irrespective of 
A. 

The galaxy spectrum may be approximately rep- 
resented as a redshifted template spectrum convolved 
with a certain broadening function caused by inner 
motions of stars in the galaxy along the line of sight: 
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Figure 2: Spectra of the template star and of the SO 
galaxy nucleus in the region of the doublet Nal ob- 
tained with the MPFS. The wavelength displacement 
of the galaxy spectrum is due to the Doppler shift, the 
line broadening is caused by the velocity dispersion 
along the line of sight. 

the broadening function B(n) in a first approximation 
can be represented by a Gaussian: 
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where a is the numerical coefficient, v and a are the 
velocity and velocity dispersion along the line of sight. 

The methods described below are based on calcu- 
lations of a discrete Fourier transform defined for the 
function F{n) as 
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where N is the number of points in the spectrum. 



2.2. Cross-correlation method 

The TD method is based on computation of a nor- 
malized cross-correlation function determined as 



C{n) = G{n)®T{n) = 

EG(m)T(ni — n), 
\ / \ 
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where 



n n 

After Fourier transformation, relationship (||) is 
written as 
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G[n) w aT{n) B{n - v), 



(3) where stars marks the complex conjugation. 



Applying the procedure of cross-correlation to ex- 
pression one can readily obtain 



C{n) = (r(n) (g> T{n)) B{n - w), 



(8) 



i.e. the cross-correlation of the galaxy and star spec- 
tra is a convolution of the auto-correlation of the tem- 
plate star spectrum T{n) (which provides information 
about the instrumental profile of the spectrograph) 
with the broadening function B{n). The velocity v is 
thus defined from the location of the cross-correlation 
function C(n) maximum. 

If the central peak of the cross-correlation func- 
tion may be represented by a Gaussian with a disper- 
sion /i: 



C {n) — Co exp 
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(9) 



and the peak of the auto-correlation function of the 
template spectra by a Gaussian with a dispersion 
^/2T■. 



T{n) (g) T{n) = exp 



(10) 



then the square of the velocity dispersion is computed 
as 

<j^^I?^2t\ (11) 

In Fig.^ are presented examples of the star and 
galaxy spectra prepared for cross-correlation, as well 
as the auto-correlation and cross-correlation peaks. 
In practice, due to the difference between the shape 
of the peak and the Gaussian one, to the contribution 
of the low-frequency component etc., the relation be- 
tween a and /i may differ from ( pi] ) , which introduces 
systematic errors into the estimate of cr. 

A number of authors have proposed different 
modifications of the TD method for estimating the 
velocity dispersion. So, in the paper by Bottema 
(1988), straightforward use of relation (^) is made — 
the template auto-correlation function is convolved 
with Gaussians of different width until the cross- 
correlation peak shape is best reproduced. We dwelled 
on the application of the empirical relation between 
the cross-correlation peak width and the velocity dis- 
persion (see Nelson and Whittle, 1995). For this the 
template spectrum T is convolved with Gaussians 
of different width ain , a cross-correlation function of 
broadened spectrum with T is constructed, and using 
formula (|l^) the velocity dispersion Oout is formally 
determined. The relationship between Uin and Uout, 
presented in Fig.|4[ is approximated by a second-order 
polynomial: 



i=0 



(12) 



and is next used to estimate the dispersion from the 
cross-correlation peak width: 

a = ya,{^?-2r^y'\ (13) 



2.3. The Fourier quotient method 

This method was put forward in the SSBS paper. By 
applying Fourier transform (^ to (H) with provision 
for (0), obtain 
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the real part in the exponent is defined by the veloc- 
ity dispersion, the imaginary part — by the galaxy 
velocity. From relation ([ij) both the parameters (a 
and v) are formally determined, however, in practice, 
the ratio of the Fourier spectra of the object and tem- 
plate is more noisy than the cross-correlation function 
peak in the TD method. For this reason, it is sim- 
pler and more reliable to find v from the shift of the 
cross-correlation peak; in this case the dispersion is 
estimated from: 
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An example of approximation of the right side of 
( |l5|) is shown in Fig.||. It is necessary to choose the 
boundaries in frequency fc, within which the fitting 
fci < A: < A:2 is made, in order to avoid the influence of 
high- and low- frequency noises. Our experience shows 
that ki = 5—10 and ^2 = 50 — 120 are optimum 
values for the number of points in the spectrum N = 
500 — 1000, being in good agreement with SSBS. 

The advantage of the Fourier quotient method 
as compared to the cross-correlation one is that to 
estimate the dispersion in the SSBS-technique, no 
additional constraints on the shape of the cross- 
correlation and auto-correlation peaks are required, 
and only the assumption that broadening function 
(^ is of Gaussian shape is applied. 

3. Instrumental profile problem 



Expression (10) shows the fact that if the spectro- 



graph instrumental profile is described by a Gaussian 
with a dispersion r, the width of the auto-correlation 
function of the spectrum produced is then ^/2 r (pro- 
vided that the lines in the template spectrum can be 
considered to be infinitely narrow with respect to r). 
The width (and even the shape) of the instrumen- 
tal profile of the spectrograph changes over the field 
(along the spectrograph slit and over the multipupil 
array), which is caused first of all by the spectrograph 
optics, precision of its adjustment and focusing. Be- 
sides, due to temperature fluctuations and flexures of 
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Figure 3: Cross-correlation method, a, b are the spectra of a Kill star and the galaxy NGC 2950, the continuum 
subtracted; c, d are the spectra on a logarithmic scale with a step of 50 km/s, the boundaries of the spectra are 
multiplied by a sinusoid (see 4-1); e, f are the auto- correlation and the cross-correlation functions, the dotted 
line shows the approximation by a Gaussian. 



the spectrograph, the quantity r at each point is a 
function of time of observations and zenith distance. 

Analysis of spectra taken with the new version of 
the MPFS shows that the spectrograph instrumental 
profile is time-stable, i.e. the variations of r during 
the observing night can be neglected. In this case, 
the profile width is defined only by the coordinates of 
the spatial element a given spectrum corresponds to: 
r = T{x,y). 

Fig.^ presents estimates of the MPFS profile 
width for different spectra (a He-Ne-Ar lamp com- 
parison spectrum, field-of-view defocused template 
star spectrum, twilight sky spectrum). The measure- 
ments made from the comparison spectrum are sys- 



tematically larger by 20-30 km/s (about 30% of the 
mean). This is connected with the fact that the en- 
trance pupil image for the calibration beam is dis- 
placed with respect to the pupil position when ob- 
serving the object. It can be seen from Fig. ^ that the 
profile width variations are ±(10-;- 20)% of the mean 
valuej] (r « 70 km/s). 

We performed experiments using different meth- 
ods of taking account of the systematic variations 
of the instrumental profile with the aim of reliable 



^ Note that all the measurements of the instrumental profile 
are made in terms of dispersion of the Gaussian. The deter- 
mination of the profile width as full width on half maximum 
FWHM fa 2.355r is more usual. 
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Figure 4: Relationship (12): am is the width of the 
smoothing Gaussian, Uout is the dispersion estimated 
from the cross- correlation peak width (dots). The solid 
line is the approximation by a second-order poly- 
nomial, the dotted line corresponds to the relation 



mapping of the velocity dispersion distribution a and 
came to a conclusion that the following procedure 
is the most correct and simple at the same time. A 
bright (5™ — 6™) template star is defocused over the 
whole array of micropupils (16"x 15"). During moder- 
ate exposures (90-200 s) at the MPFS it is possible to 
take template spectra with high signal-to-noise ratio 
over the whole field. When analysing the galaxy spec- 
trum by the TD and SSBS at each point of (x, y), an 
individual template spectrum T{x, y) is taken. Such 
a procedure takes actually full account of the effect of 
the instrumental profile, provided, of course, that the 
spectrograph is carefully focused before the observa- 
tion. 



4. Sequence of work with MPFS data 
4.1. Preparation of spectra 

Using the language IDL, we have written a package of 
programmes CROSSLIB which implements the meth- 
ods presented in Sections 2 and 3. The initial data 
are 2D spectra of a galaxy and a star (Fig. |l]) after 
preliminary reduction (bias correction, flat-fielding, 
cosmic hits removal, individual spectrum extraction 
from the CCD images and its conversion to a com- 
mon linear wavelength scale, night-sky line subtrac- 
tion, correction, if needed, for the effect of differential 
atmospheric refraction) . 




Figure 5: Approximation of function (15) by a Gaus- 
sian in the Fourier quotient method. The boundaries 
of fitting (ki and k2) are marked. 



After the preliminary reduction the continuum 
is subtracted from the spectra. Underestimation of 
the continuum level leads to appearance of a low- 
frequency component in the Fourier spectrum and, 
as a consequence, to greater noises for small k in 
the SSBS method, and of a wide base in the auto- 
correlation and cross-correlation peaks in the TD 
method. It is not infrequent that high-order polyno- 
mials are used to fit a continuum (Larsen et al., 1983). 
However, our experience of reduction has shown 
that the best result can be achieved if a spectrum 
smoothed by a median filter with a window of 150- 
200 pixels is taken as a continuum. In so doing, even 
the broadest absorption spectral features are practi- 
cally undistorted. 

The emission lines in the galaxy spectrum reduce 
the contrast of the cross-correlation peak. Depending 
on the number and the intensity of lines, different 
techniques can be employed to remove them. In most 
cases it suffices to approximate by a Gaussian the 
brightest line, determine the intensity of the rest of 
the lines with respect to it and subtract from the 
spectrum the Gaussians with corresponding weights. 

After the continuum subtraction from the emis- 
sion lines, the spectra are rebinned onto a logarith- 
mic scale according to (|l]). To diminish the effect of 
"leakage of energy" in the Fourier transform due to 
the sharp cutoff of the ends of the spectrum (for more 
details see Brault and White, 1971; TD, 1979), the 
spectrum needs to be multiplied by the function of 
the form: 
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Figure 6: Instrumental profile of the MPFS (dispersion of the Gaussian t in km/s) evaluated in the spectral 
range 5100-6100 A, the spectrograph dispersion is 1.35 A/px: a — from the comparison spectrum, b — from 
the spectrum of the Kill star, c — from the twilight spectrum. 



here p = 0.05 — 0.1 is the relative size of the spectrum 
boundaries multiplied by a sinusoid; examples of such 
spectra are displayed in Fig.||c,d. 

4.2. Determination of v and a 

For each spectrum corresponding to one of the spatial 
elements {x,y), relation ( |T^ ) is constructed. For this 
purpose, the spectrum of the template star T{x, y) is 
smoothed by Gaussians of different widths and an es- 
timate of the dispersion is made by the TD method. 
Because the cross-correlation peak shape is not al- 
ways precisely Gaussian, it is necessary to fix the 
boundaries of the region of fitting ( [uq — An, Uq + An] , 
where Uo is the peak location), in which a robust fit 
with a Gaussian is executed, considering that calibra- 
tion relation (^) may change for different An. 

Obtain a set of coefficients: ^{x,y), ai{x,y), i — 
0,1,2. If the spread of the coefhcients from point 
to point is considerable, it is better to approxi- 
mate these coefficients along the slit by a second- 
order polynomial. Construct cross-correlation func- 
tions G{x,y) '^T{x,y) using relation (^), determine 
V, and from the peak width determine a with the 
from of (|ll|) and (|l|) (Fig.|). One should use here 
the same value of An as the one used in construction 
of calibration relation ( p^ ) . 

To estimate the dispersion by the SSBS method, 
the right side of expression ( p^ ) is approximated with 
allowance made for the remarks on the choice of the 
limits fci and k2 (see subsection 2.3 and Fig. ||). 

By applying the above-described procedure to 
the spectra of all spatial elements, produce two- 
dimensional maps of radial velocities v{x,y) and ve- 
locity dispersions aTD(x,y), (yssBs{x,y) (the index 
indicates the method of construction). The heliocen- 
tric velocities (reduced to the centre of the Sun) are 



found from the expression: 

VHel=V + V,-VQ,,+VQ^gal. (16) 

Here F©,*, V^^gai are the heliocentric velocity correc- 
tions for the template star and the galaxy computed 
by the standard method, K is the radial velocity of a 
star from the catalogue (reduced to the centre of the 
Sun) . 

4.3. Choice of spectral range 

The optimum spectral range, in which the analy- 
sis described above is made, must contain absorp- 
tion lines characteristic of the stellar population of 
galactic nuclei, unblended, when possible. It is de- 
sirable that the number of emission lines of ionized 
gas should be a minimum. Most commonly the re- 
gions around the triplets Mgl AA 5167, 5172, 5183 A 
and Call AA 8498, 8542, 8662 A are used. The latter 
is more preferable because with one and the same 
dispersion of the spectrograph the spectral scale (in 
km/s) in the region of Call is 1.6 times as large as in 
the region of Mgl. For the advantages of using Call 
see the paper by Terlevich et al. (1990). 

Unfortunately, the available set of gratings in 
MPFS made it impossible to work in the region of 
the infrared triplet Call without substantial losses 
of light. Experiments performed in the optical range 
have demonstrated that the most optimum spec- 
tral range to study kinematics of stars at the cen- 
tres of SO~Sb galaxies is 5100 - 6100 A. It includes 
strong absorption hues of Mgl A 5175 A, Fel A 5229 A, 
FeI-hCaIA5270A, NaIA5893A and others. At the 
same time, the emission lines are relatively weak and 
not numerous, noticeable mainly only in Sy galax- 
ies ([NI] A 5199 A, [FeVII]AA 5158, 5276, 5720, 6086 A, 
Hel A 5875 A) . The observations are made with a grat- 
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Figure 7: Errors of determination of v ( top ) and a (in the middle by the SSBS method, at the bottom by the 
TD method) for the galaxy NGC2950 over the whole MP FS field of view. The errors were estimated from 4 
template stars. The left-hand panels are the histogram of distribution of errors and their Gaussian fit. The 
right-hand panels show the errors as a function of cross- correlation peak height. 



ing of 1200gr/mm yielding a reciprocal dispersion of 
1.35A/px. The spectral range 4800-6100A contains 
both the absorption lines described above and the 
brightest emission lines H^, [OIII] A 4959, 5007, which 
allows a simultaneous study of ionized gas kinematics. 

4.4. Assessment of errors 

The assessment of errors in different methods of de- 
termination of (T and V is considered in some detail in 
a number of papers (TD, 1979; SSBS, 1977; Larsen 
et al., 1983). The influence of random and system- 
atic errors is, however, best seen when comparing the 
results obtained from different template stars during 
one night. Fig.^ shows the distribution of errors in 



the determination of the velocity and velocity disper- 
sion for all 240 individual spectra of NGC 2905 in the 
MPFS field of view. The mean value for each spa- 
tial element was estimated from four stars of spec- 
tral classes G7 III-K2 III. The values of the mean er- 
rors estimated by the two methods, as the second 
moment in the distribution of deviations of individ- 
ual measurements from the mean value over the four 
stars and as the width of the Gaussian that approxi- 
mates this distribution, are indicated beside the his- 
tograms of the errors. These errors are typical of our 
observations and make Sv — 5 — lOkm/s for the 
velocity, SaxD = 10 — 15km/s for the velocity dis- 
persion estimated by the cross-correlation technique. 
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Figure 8: Comparison of systemic velocities (a) 
and central velocity dispersions (b) derived from the 
MPFS observations with the data of the catalogue 
RC3 and the HYPERCAT database. 

and SassBS = 15 — 25km/s for the Fourier quotient 
method. For the central regions where the value of the 
cross-correlation peak is Cmax = 0.8 — 1.0, these er- 
rors are about 1.5 times as small as the field average. 
It can be seen from Fig.0 that the velocity disper- 
sion errors in the SSBS method increase sharply at 
Cmax < 0.7, that is, one can estimate reliably only 
the central velocity dispersion. However, since this 
method requires no additional calibration (subsection 
2.3), it is good practice to use the SSBS method esti- 
mates for checking the absolute values of the central 
velocity dispersion. 

We have carried out a number of experiments 
for testing the precision of the TD and SSBS tech- 
niques, using for this purpose the spectra of stars 
of class G6III/IV-K3III obtained with the MPFS. 
These spectra were artificially smoothed and ren- 
dered noisy, then they were employed to estimate the 
velocity dispersion values. Analysis has shown that 
if the signal-to-noise ratio in individual spectra is 
S/N > 20, it does not affect the accuracy of a deter- 



mination, which is less than 5 — 10% of the absolute 
value of ct; this is in good agreement with Fig.^. 

The experiments with smoothing have confirmed 
the well-known fact that it is impossible to reliably 
determine the velocity dispersion if it is smaller or 
close to the instrumental profile width. In our case 
this limit is 70-80 km/s (see Section 3). 

The main factor that affects the accuracy of de- 
termination of the dispersion for a given S/N and 
spectral interval are the differences in spectral fea- 
tures of a galaxy and a star which lead to decreasing 
the contrast of the cross-correlation peak and to dif- 
ferences in the Fourier spectra of the SSBS method. 
When studying the motions of stars in galaxies of 
early types, the spectra of red giants G6III-K3III 
are normally used as templates. It is a good idea to 
take spectra of stars of different spectral classes from 
this range during one observing night to select a tem- 
plate which is the most close to the galaxy spectrum. 

For 8 galaxies observed with the MPFS, we com- 
pared our measurements of the systemic velocity 
(the galaxy centre velocity) with the data of obser- 
vation of neutral hydrogen velocities in the 21cm 
line, which had been taken from the catalogue RC 3 
(Fig.||a). The root-mean-square value of deviations 
was 17km/s without a systematic shift. Data on the 
nuclear velocity dispersion of stars were found in the 
database HYPERCAT]^ and compared with our es- 
timates (Fig.^D). No systematic differences are ob- 
served; the great scatter of points is most likely due 
to the fact that only using the 2D technique, one can 
reliably subtract the spectra from the galaxy dynam- 
ical centre. 

5. Example of data reduction: 
NGC2950 

According to Wozniak et al. (1995) and Friedly et 
al. (1996) the SBO galaxy NGC2950 refers to the so- 
called double-barred galaxies. In the optical and IR 
images of the galaxy there is a bar of about 6" in 
size inside a 38" primary bar. The position angles 
of the bars are PA = 332° and 272°, respectively. 
NGC2950 was observed with the 6 m telescope on 
March 27/28, 2000 using the spectrograph MPFS, 
the spectral range was 4820 — 6200 A, the dispersion 
1.35 A/px, the scale 1" per lens, the seeing 2". The to- 
tal exposure was 3600 s. The field of view of 16" x 15" 
covered a region of the inner bar and the central part 
of the outer bar. 

Prior to the continuum subtraction, a smoothing 
of data was performed in each spectral bin by a Gaus- 
sian of FWHM — 1" X 1" with the aim of optimum 
filtering of noises. In so doing the spatial resolution 
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Figure 9: MPFS maps for the galaxy NGC 2950. a — the continuum image (the surface brightness logarithm 
is in ergs/crn^ /s/D" ), the stright line shows the orientation of the main (outer) bar; b — the velocity field of 
stars; c — the distribution of the velocity dispersion constructed by the TD method. The origin of coordinates 
is at the dynamic centre of the galaxy (shown by cross). 



decreased slightly and was about 2'.'3. The further 
processing was conducted in accordance with Section 
4. After the rebinning onto a logarithmic scale, the 
velocity bin was 50 km/s. Five template stars were 
taken during the night, 4 of which, yielding a maxi- 
mum correlation with the galaxy spectrum, were used 
to construct average velocity fields and velocity dis- 
persion (Fig.||). When constructing these fields, the 
points where the error was larger than 20 km/s were 
masked. 

The central velocity dispersion values measured 
by the TD method {a = 205 ± 5 km/s) and by 
the SSBS method (a = 190 ± 10 km/s) agree with 
each other and with the data of HYPERCAT (170 ± 
10 km/s, 182 ±14 km/s and 210± 14 km/s from differ- 
ent estimates). The central ellipsoidal structure elon- 
gated along PA 330°, i.e. the one associated with 
the "outer" bar, is distinguished in the distribution 
of velocity dispersion (Fig.^). At the present time, 
we continue the work over the study of kinematics of 
double-barred galaxies; preliminary results are avail- 
able in the paper by Moiseev and Afanasiev (2001). 

6. Conclusions 

With the above-described procedure one can con- 
struct velocity fields and velocity dispersions of stars 
in the central parts of early-type galaxies on the ba- 
sis of data obtained with the spectrograph MPFS. 
The accuracy of velocity measurement is on the av- 
erage 6v = 5 ~ 10 km/s, that of velocity dispersion 
measurement is Saro = 10 — 15 km/s and SassBS = 
15 — 25 km/s for TD and SSBS methods, respectively. 
This technique allows us to take account of the varia- 
tions of the instrumental profile and produce reliably 
two-dimensional maps of the velocity dispersion if it 



is not less than 70-80 km/s. The systemic velocities 
and the central velocity dispersions that we have mea- 
sured in 8 SO~Sb galaxies and the HYPERCAT data 
do not show systematic differences. 

We associate the future prospects of studying 
kinematics of stars in galaxies aided by the spec- 
trograph MPFS with the possibility of observing in 
the region of the infrared triplet Call. The accuracy 
of measurement of velocities and velocity dispersion 
must in this case be about 1.5-2 times better. 
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